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Abstract. We study the evolution of scalar curvature perturbations in a brane- 
world inflation model in a 5D Anti-de Sitter spacetime. The inflaton perturbations 
are confined to a 4D brane but they are coupled to the 5D bulk metric perturbations. 
We numerically solve full coupled equations for the inflaton perturbations and the 5D 
metric perturbations using Hawkins-Lidsey inflationary model. At an initial time, we 
assume that the bulk is unperturbed. We find that the inflaton perturbations at high 
energies are strongly coupled to the bulk metric perturbations even on subhorizon 
scales, leading to the suppression of the amplitude of the comoving curvature 
perturbations at a horizon crossing. This indicates that the linear perturbations of the 
inflaton field does not obey the usual 4D Klein-Gordon equation due to the coupling 
to 5D gravitational field on small scales and it is required to quantise the coupled 
brane-bulk system in a consistent way in order to calculate the spectrum of the scalar 
perturbations in a brane-world inflation. 
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1. Introduction 

Our view of how our Universe may be described within a higher-dimensional spacetime 
has been revolutionised via brane-world models [1] (see, e.g., [21 [3j H] for review). A 
simple way to study the effects of the higher- dimensional gravity on inflation models 
is to consider a 4D inflaton field confined to the brane in the context of the Randall- 
Sundrum model [5j. Work to date has concentrated mainly on deriving corrections to 
the standard 4D result by assuming that 5D metric perturbations generated by the 
inflaton perturbations can be neglected. Then the modification of spectrum is coming 
from the modified Friedmann equation. Although a lot of interesting results have been 
obtained in this direction [6j [7], it remains to see the consistency of neglecting the bulk 
metric perturbations. The analysis of the evolution of gravitational waves in this model 
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has revealed that the creation of massive graviton modes in the bulk is as important as 
the modification of the Friedmann equation [Sj. In fact, these two higher-dimensional 
effects cancel out to give a nearly identical spectrum with the conventional 4D models 
for the stochastic background of gravitational waves in the radiation dominated universe 
[9]. This highlights the importance of taking into account the higher-dimensional effects 
consistently. 

The study of the effects of bulk metric perturbations has been initiated in 
Ref. pHl E] (see [12] for a different approach). Using the approximation that the 
geometry of the brane is approximated by de Sitter spacetime for a slow-roll inflation, the 
bulk metric perturbations sourced by inflaton fluctuations on the brane were calculated. 
It turns out that the inflaton fluctuations excite an infinite ladder of massive modes of 
bulk metric perturbations. Then Ref. [TT] tried to estimate the effects of the backreaction 
of bulk metric perturbations. Assuming that the inflaton field behaves in the same ways 
as 4D models on sufficiently small scales at the 0-th order of slow-roll parameters, 
the corrections to the evolution of inflaton perturbations at the first order of slow-roll 
parameters were calculated. Unlike in 4D models where the slow-roll corrections can 
be completely ignored on sub-horizon scales, the slow-roll corrections play a role even 
on sub-horizon scales, yielding the change of the amplitude at a horizon crossing. This 
result may not be so surprising as the bulk metric perturbations are strongly coupled to 
inflaton perturbations on small scales where gravity becomes 5D. However, this indicates 
that we should carefully check the validity of using the 4D formula just by modifying 
the background Friedmann equations. 

This paper focuses on the classical evolution of inflaton perturbations coupled to 
bulk metric perturbations. We do not assume de Sitter geometry on the brane and 
we take into account the backreaction of inflaton dynamics consistently. This entails 
a numerical analysis to solve the coupled system directly. We investigate whether the 
inflaton perturbations behave as free massless fields on small scales. 

The structure of this paper is as follows. In section II, we set up basic equations 
to describe the evolution of inflaton perturbations coupled to metric perturbations. 
We use the Hawkins-Lidsey model of the inflation on the brane where the background 
solution is derived analytically [13]. The Gaussian-normal coordinate is used to solve 
the perturbations. Section III is devoted to the numerical scheme. In section IV, we 
explain initial conditions and boundary conditions. Then we study the evolution of the 
comoving curvature perturbations, which are directly related to observables. 

2. Basic equations 

2.1. Background spacetime 

We consider the Randall- Sundrum single-brane model [14]. In this model, a 3D spatial 
hypersurface is embedded in a 5D anti-de Sitter (AdSs) space-time with a negative 
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cosmological constant A 5 . The 5D Einstein equations in the bulk are given by 

^G AB + A 5 (5 W = 0. (1) 

The cosmological constant is related with the curvature scale of the AdSs space-time 
\x as A 5 = — 6/x 2 . The embedded brane has a tension A given by k 2 X = 6/i. Here k 5 
denotes the 5D gravitational constant which can be rewritten using the 5D Planck mass 
M 5 as k\ = 87r/Mf. On the other hand, we use K4 as 4D one defined as n\ = Sn/M^ 
which is related with k 5 as k 2 = finf, where M p \ is the 4D Planck mass. The induced 
metric on the brane is defined as 

9ab = {5) 9ab ~ n A n B , (2) 

where n A denotes the unit normal vector to the brane. The energy-momentum tensor 
T^ u on the brane and a brane tension —\g^ u produce the jump in the extrinsic curvature 
at the brane; 



T\ - ljf v {T - A) 



(3) 



where we imposed the Z2 symmetry on the brane and the extrinsic curvature is defined 
by K^ v = g^ c g u D ^V c n D- According to Ref. [15], the effective Einstein equation for 
the induced metric on the brane is obtained as 

G/j,u = K^r^ v + K^n^y e^v, (4) 

where IT^ is the quadratic terms of the energy-momentum tensor, 

111 1 

IV = ~^ T ^T u a + y^TT^ u + -g^ u T af3 T a P - —g^T 2 , (5) 

and is the electric part of 5D Weyl tensor Cabcd, 

E^ v = C E ' AFBn E n F g^ A g u B '■ (6) 

We assume that the matter content confined to the brane is described by a perfect 
fluid whose energy-momentum tensor is given by 

r"„ = diag(-p,p,p,p), (7) 

where p and p denote the energy density and the pressure of matter on the brane, 
respectively. In this paper, we use the Gaussian-Normal (GN) coordinate. The metric 
is given by pa] 

ds 2 = -n 2 (y, t)dt 2 + a 2 (y, t)5 lj dx i dx j + b 2 {y } t)dy 2 . (8) 
Here the warp function a(y,t) and the lapse function n(y,t) are obtained in [16] as 

a(y, t) = a (t) jcosh fiy - ^1 + sinh fiy J , (9) 
n(y,t) = cosh jiy - |l - ^(2p + 3p)| sinh/!?/, (10) 
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where a$ is the scale factor on the brane, and we can set b(y, t) = 1 without a loss of 
generality. The scale factor is determined by the Friedmann equation and the energy 
conservation law as 

K l„ , P 



P =S(p + p)H. (12) 

In this coordinate, the brane is located at y — — 0. Note that, in this coordinate, 
there is a coordinate singularity at the finite distance from the brane, 

y c = - coth- 1 ( — ) , (13) 
P \P J 

where the warp function vanishes a(y c , t) = 0. The location of the coordinate singularity 
y c agrees with the Cauchy horizon of the AdSs spacetime [T7]. We will come back to 
this issue later. 

In this paper, we consider a scalar field on the brane. The scalar field satisfies the 
equation of motion 

4> + 3H<p + V'(<f>) =0, (14) 

where V(<f>) denotes the potential of the scalar field. We use a model proposed by 
Hawkins-Lidsey as an inflaton model on the brane to avoid numerical integrations of 
the background equations p3]. The potential for inflaton is given by 

V(4>) = ^(6 - C 2 )cosech 2 {^-^ , (15) 

where C is constant. An advantage of this potential is that an exact solution for the 
scale factor is known 

a (t) = [(CV + lf-l] 1/C2 . (16) 

The energy density and pressure of the scalar field can be calculated as 

2A 

P (CV+l) 2 -!' 

O \ 1 

^T' c2 - 3 > (CV + D--r (18) 

Thus the effective equation of state w becomes 
r> r 2 

w = P - = -l + ( -. (19) 
P 3 

Then for small C 2 , a slow-roll inflation is realized, which is a generalisation of the 
power-law inflation. 
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2.2. Scalar perturbations 

2.2.1. Equations in the bulk Now let us consider the scalar perturbations. There are 
several ways to calculate the scalar perturbations in a AdSs spacetime. Here we describe 
an approach based on a master variable [111 [19]. Using the generalized 5D longitudinal 
gauge, the perturbed spacetime is given by 

ds 2 = - n(y, t) 2 (l + 2A)dt 2 + a(y, t) 2 (l + 2K)5 ij dx i dx j 

+ (1 + 2A yy )dy 2 + n(y, t)A y dydt. (20) 

It was shown that the solution for metric perturbations can be derived from a master 
variable Q as 

2Q!' - -n') + f n - ) - rfn } , (21) 

n' 
n 



A 


1 

6a 


Ay 


= ±{ 


na \ 


A 


6a [ 


K 


6a I 



n l n 2 \ n 



n',\ 2 I - n ■ 



(22) 

— \ [q!'-2—q! ) + 4 [n-'-n] +i2 2 n\ , (23) 

' \ n J n 2 \ n J J 

Uf + Vlq>\ + L f-n + ^ _ 2/ /o j , (24) 



as long as the master variable Q in the bulk satisfies a wave equation given by 
1 aY / n _,\ ' / o V 2 \ n 



na 3 J \o? ) \ a 2 J a 3 



n + K n ' +U 2 + — -fi = o, (25) 



where a prime and a dot denote the derivatives with respect to y and £, respectively. 

2.2.2. Boundary conditions In this section, we derive the evolution equation for 
inflaton perturbations coupled to the master variable. See Refs. [TQI EDI ED E2[ E3] 
and references therein for earlier works. We begin with the effective Einstein equations 
on the brane [IS]. The perturbed effective Einstein equations are given by 

(4) 8G„ V = KlSTp, + -5E lMU . (26) 

The perturbed Einstein tensor is given by [21] 

2 

al 

^8G\ =-2(HA h -TZ h ) tt , (2* 



^SG° = QH (-K h + HAC) + ^V 2 ^ b , (27) 



(29) 



where 



{4) 5G T = 2 { (3H 2 + 2H)A h + HA h -H h - 3HU h + -L V 2 {TZ h + A h ) \ , (30) 
I 3a^ J 

{4) 5G TF = -\(n h + A h ). (31) 



a 2 
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The subscript b denotes a quantity evaluated on the brane. We can parameterize the 
projected Weyl tensor 5E^ U as an effective fluid; 

5E% = KlSpe , (32) 

5E\ = -Kl8q £ji , (33) 

6E l 3 =-K\{8p e 8r i j + toej) , (34) 

where 

Sir £ l j = Sirs' 1 j - -^S l j8irs' k ,k ■ (35) 

The Weyl fluid satisfies the radiation-like equation-of-state 5ps = Sps/3. The perturbed 
energy-momentum tensor ST^ is defined as 

ST° = -Sp , (36) 
ST° t = Sq A , (37) 
STj = Sp&j . (38) 

Throughout this paper, we focus on the perturbations of the scalar field. Thus the 
anisotropic stress 57r l j on the brane vanishes. The quadratic energy-momentum tensor 
SH^ is calculated as [21] 

511% = -\p8p , (39) 
6 

SU°i = lpS q>i , (40) 
o 

Sn*i = ~{(p + p)Sp + P Sp}S i j ■ (41) 

All perturbative quantities can be expanded by scalar harmonic functions e lkx , which 
enables us to replace the operator V 2 by —k 2 . Then, from the effective Einstein 
equations (126]) . we can derive the evolution equations for the perturbative quantities 
as [23] 

1 k 2 1 

- -4mp = 3H{n h - HA h ) + —n h - -4s P s, 

—KfHSp =-ti h -H U-IL\n h + HA h + [2H + 3H 2 



-KpiSq 




b 




a 2 

where we defined 
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and used the relation p + p = 2H/k\H. derived from Eqs. (TTTT) and (TT2"|) . Ref. [22J showed 
that the components of the perturbative Weyl tensor 5p£,5q£,5iT£ can be written in 
terms of the master variable Q as 
2 r ( k 



K\8q £ 



3^/ b 
k 2 

~ ( HVt - 

3oq V / 1, 



6oq [ Oq \ n a 



(47) 
(48) 
(49) 

(50) 
(51) 
(52) 

The equation of motion for the scalar field perturbations is derived from the conservation 
as 



The perturbed energy-momentum tensor is calculated as 
6p = -<p 2 A h + 050 + V'((f>)S4>, 
6p = -<p 2 A h + j>6<j> - V'((f>)6<f), 
8q = -<; 



law V^T^ 



k 2 

5(p + 3H5(p + V 



-V"(<t>)6<t> - 30ft b + <j>A h - 2V'(<f>)A h . 



(53) 



Now we introduce a gauge-invariant variable, the Mukhanov- Sasaki variable, defined as 



Q 



H 



(54) 



Substituting this definition into Eq. fl5"3l and using the effective Einstein equations, we 
can derive the evolution equation for Q as 

H IIY'(o) J hY , t „. 



Q + 3HQ + ^Q+r--2 

a I si si 



-2(-| +V"{<P)}Q = J(fi), 



(55) 



where 



.(56) 



Here the left-hand side of Eq. f !55p is completely the same as the 4D cosmology while 
J(Q) describes the the effect of the bulk metric perturbations. Using Eq. (T251) . the 
curvature perturbations ( 1241) evaluated at the brane 71^ can be expressed in terms of Q 
and its first-order derivatives as 



If k 2 
3HQ' - 3HQ - 3fi 2 Q + — Q 



(57) 



6ao V a o 
From the (0, i) component of the effective Einstein equation fj4~4j) . we can relate Q and 
Q as 

H 



Q 



Q' -HQ 



HH 



{ k 2 

nn' - u 2 n + -^n 

V 3a o 



(58) 
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which gives a non-local boundary condition for Q. Alternatively, we can combine 
Eqs.(PJ and (jHJ to get 

a 4(5p-3H5q) = (n' --n) . (59) 

Then using Eqs.( l50p and fl52|) and rewriting 8(f) by Q, we obtain 

K ^ 2 (f Q ) ' = "5 {ir (« - « n ) +H - } b • < 60 > 

In summary, the master equation (125]) . the Mukhanov- Sasaki equation ([55]) with a source 
( 156]) and the junction condition (|58]) or ( 160]) form a closed system for the coupled inflaton 
and metric perturbations. 

Note that, combining Eqs. (142]) - (jHJ), one can derive a boundary condition for Q 
which does not include Q [25j . This boundary condition is equivalent to choose ( 155]) 
and (158]) as boundary conditions. 

3. Numerical scheme 

In this subsection, we briefly discuss the numerical scheme to solve the equations. First 
of all, in order to avoid the effect from the coordinate singularity at y = y c in the 
GN coordinate, we introduce an artificial cutoff (regulator) boundary in the bulk at 
V — Dreg = IVc- The constant 7 controls the position of the regulator brane. Since the 
regulator brane should be far from the physical brane and 7 < 1, we set 7 = 0.99 in 
our numerical simulations. We discuss the boundary condition at the regulator brane 
in the next section. 

In this paper, the numerical calculations of this system were carried out by 
employing the pseudo-spectral method [26]. The master variable Q(t,y) is decomposed 
by the Tchebychev polynomials defined as 

T n (0 = costncos^O for - 1 < f < 1, (61) 

which is a polynomial function of £ of order n. Here the variable £ is related to the GN 
coordinate y. To implement the pseudo-spectral method, we must transform a spatial 
coordinate y in the GN coordinate by 

y = g { fa** -Vb)Z+ (yveg + yb)}, (62) 

so that the locations of both the physical and the regulator branes are kept fixed and 
the spatial coordinate y is projected to the compact domain — 1 < £ < 1. Adopting this 
coordinate system, the quantity Q(t, £) is first transformed into the Tchebychev space 
through the relation, 

N 

fi(t,£) = ^fi n (t)T n (£), (63) 

n=0 

where we set iV = 64 or 128. We then discretize the £-axis to the iV+l points (collocation 
points) using the inhomogeneous grid £ n = cos(n7r/iV) called Gauss-Lobatto collocation 



Numerical study of curvature perturbations in a brane-world inflation at high- energies^ 



points. With this grid, the fast Fourier transformation can be applied to perform 
the transformation between the amplitude Q(i, £) and the coefficients Q n (t). Then 
the master equation ( 1251) rewritten in the new coordinates are decomposed into a set 
of ordinary differential equations for Q n (t). For the temporal evolution of and 
Q{t) satisfying Eq. (l55l) . we use Adams-Bashforth-Moulton method with the predictor- 
corrector scheme. Further technical details of the numerical scheme is summarized in 
the appendix of Ref. [9]. 

4. Evolution of curvature perturbations 

4-1. Initial conditions and boundary conditions 

The boundary condition at the brane is given by the junction condition Eq. ( 1581) or 
Eq. (160]) . At the regulator brane, we use a Dirichlet boundary condition 

^U rcg = 0. (64) 

This is because the analysis in de Sitter brane background pJJ, [UJ suggests that the 
inflaton fluctuation excites the bound states that decays quickly away from the brane. 
If we focus on this bound state the result should not be sensitive to the location of the 
regulator brane. We will check this fact numerically later. Different boundary conditions 
from Eq. (I64|) may produce reflection waves from a regulator brane, which affect the 
evolution of Q. A proper boundary condition should be imposed on the past Cauchy 
horizon. This requires a quantisation of 5D metric perturbations in the bulk, which 
will turn out to be a highly non-trivial problem due to the coupling to the inflaton 
perturbations. 

In 4D cases, Q satisfies 

k 2 

Q + 3HQ + —Q = 0, (65) 
a 

on small scales since k 2 term is dominant over the mass term induced by the dynamics 
of the inflaton [see (1551) ]. Thus the solution on small scales is given by simple plane 
waves and we can choose the usual 4D Bunch-Davis vacuum state, which specifies an 
appropriate phase and normalization. On the other hand, in the present 5D case, it 
is not guaranteed that the perturbed inflaton field obeys the usual 4D Klein-Gordon 
equation due to the coupling to the bulk metric perturbations, which is described by 
J(Q) in the right-hand side of Eq. fl55|) . In order to obtain the appropriate initial 
conditions for Q(t), we again need a quantisation of the coupled system of the inflaton 
perturbations on a brane and 5D metric perturbations in the bulk. 

The aim of this paper is to check whether we can neglect J(Q) on small scales or 
not. If this term could be neglected on the small scales, the quantity a Q would behave 
just as plane waves with a constant amplitude as in the standard 4D cases. We call 
the solution of Eq. ( 1551) with J(f2) = or, equivalently, Eq. ( 1651) on small scales 4D 
solutions. If the solutions for the coupled equations agree with 4D solutions, we can 
assume that the inflaton perturbation behaves as a free massless field without taking 
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into account the effect of 5D metric perturbations, which is the assumption made in 
literatures. 

For this purpose, we take the simplest possible initial conditions for Q(y, t) and then 
specify initial conditions for Q(t) so that they are consistent with the junction condition 
(IBDj) . We set Q(y,ti) = 0, Cl(y, tj) = and Q(U) = 1. Then the time-derivative of Q(t) 
is determined from the junction condition fl60l) . 

Note that our initial conditions do not give standard qunatum vacuum state for 4D 
solutions. They are merely references to see the effects of J(Q) on a classical evolution 
of Q(t). 



4-2. Evolution of curvature perturbations 

We solve numerically the wave equation for the bulk metric perturbations fl25l) in the GN 
coordinates with the junction condition fl60|) supplemented by the evolution equation 
for inflaton perturbations ( |55l) . We set the parameter of the Hawkins-Lidsey model as 
C = 0.1 and the initial time for the simulations as \xti = 40. At the initial time, the 
Hubble scale is given by 

~ 2.92. (66) 

Thus we are considering high-energy region H/fi > 1. 

The observable is the comoving curvature perturbation defined by 

1l c = -S-Q. (67) 



We focus on the dynamics of this quantity. Fig. Q] shows the behaviour of curvature 
perturbations for kj\i = 5 which corresponds to k/a (ti)fj, ~ 296. The curvature 
perturbations 1Z C becomes constant on super-horizon scales, which has been shown to 
be valid even in brane world models [27]. On sub-horizon scales, while TZ C oc 1/ao 
in the 4D cosmology, a suppression of the amplitude is observed in the brane-world 
model, which is due to the coupling to the bulk metric perturbations. In Fig. [2j 
we compare the amplitude of 1Z C obtained in numerical simulations with the one 
obtained by neglecting the coupling to the bulk metric perturbations, i.e. J (CI) = 
in Eq. fl55l) . While the difference is very small for the long- wavelength modes (left 
panel; k/ao(ti)fj, = 296), the suppression becomes significant for the short- wavelength 
modes (right panel; k/ao(ti)fi = 2960). Due to this, the spectrum of 1Z C just before the 
horizon crossing acquires a scale dependence as is shown in Fig. [3j Perturbations with 
larger wavenumber stay on sub-horizon scales for a longer time than those with smaller 
wavenumber, so they receive more suppression. 

The suppression of the amplitude 1Z C under horizon may be understood as the 
excitation of metric perturbations. The suppression of the amplitude of Q is transferred 
into the enhancement of the metric perturbations fl in the bulk, which is shown in the 
left panel of Fig. HI 
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Figure 1. The evolution of curvature perturbations 1Z C in the inflationary epoch. We 
set C = 0.1, /iti = 40 and k/fj, = 5 which corresponds to k/ao(ti)fi w 296. 




Figure 2. The curvature perturbations (multiplied by the scale factor) in the 
inflationary epoch for a long-wavelength mode {left; fc/ao(tj)/i ~ 296) and for a short- 
wavelength mode (right; k/ao(ti)fj, w 2960). The solid lines represent numerical results 
and the dashed lines show the 4D predictions obtained by neglecting the coupling to 
the bulk metric perturbations, that is, by solving Eq. ([55]) with J(f2) = 0. We set 
C = 0.1 and fiU = 40. 

Here we comment on the reason why we evaluated 1Z C just before the horizon 
crossing. In the 4D case, we can choose the Bunch-Davis vacuum as an initial condition. 
Hence we can determine the amplitude and phase of the inflaton perturbations and 
obtain the spectrum of TZ C without any ambiguities. On the contrary, in the present 5D 
case, the initial conditions cannot be determined without quantising the coupled system. 



Numerical study of curvature perturbations in a brane-world inflation at high-energiesYl 




1f) -i I , , , I 

10 1 10 z 10 3 10 4 10 5 

k/arftjll 

Figure 3. The scale dependence of the curvature perturbations evaluated just before 
the horizon crossing. The horizontal axis represents the physical scale of perturbations 
evaluated at the initial time //i, = 40. We estimated the ratio |7?.^ D /7\!^ D | for each 
wave number where 7?.^ D is the curvature perturbation in the brane world model and 
TZf^ is the one in a 4D model with the identical background dynamics. 



Even in the classical level, the junction condition (158j) or ( 1601) restricts our choice of 
the initial conditions. On superhorizon scales, the initial phase of the perturbations 
essentially determines its final amplitude after the horizon crossing because the phase 
determines the amplitude of the growing mode solution on large scales. Then we get 
an unphysical k-dependence of the curvature perturbations evaluated on superhorizon 
scales due to the restricted choice of dQ/dt for given k. On the other hand, the initial 
phase is not important when we evaluate the amplitude of the curvature perturbations 
on small scales. Our aim is to see how J(Q) term changes the amplitude of Q in 
comparison with the 4D solution for various k on small scales. Therefore we chose to 
calculate the spectrum below the horizon scales. 

4-3. Check of numerical calculations 

In this section, we first check the assumption of introducing a regulator. As expected 
from the previous work, Q decays fast away from the brane as is shown in the right panel 
of Fig. HI This justifies the use of boundary condition Q(y Teg , t) = 0. In addition, thanks 
to this behaviour of Q, the evolution of the curvature perturbations does not depend on 
the location of the regulator. In fact, even at late times when the information from the 
regulator brane comes into the brane, the behaviour of the curvature perturbations does 
not change. Thus our result is not obscured by the regulator brane nor the existence of 
the coordinate singularity of the GN coordinate. 

Finally, we check the accuracy of our numerical simulations by checking the 
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Figure 4. Left : The bulk metric perturbation fl evaluated on the brane. Right : 
The snapshot of at fj,t = 41.0. At this time, the coordinate singularity in the bulk is 
located at ^j/ c (41.0) » 0.34. We set C = 0.1, ^ = 40, k/a Q {ti)u = 2960. 



constraint equations derived from the Einstein equations. These equations can be 
generally expressed as 

C(t) = J2Ck(t)fk(t) = 0, (68) 

k 

where 

f k (t) = {H h ,'JZ h ,'JZ h ,A h ,A h ,5p £ ,6q£,5TT £ }, (69) 

and Ck(t) is a set of their coefficients. We calculated the relative error for the constraint 
equations (HI) flHD and (J45]) as 

*W = \rl\ f mi • ( 7 °) 

maxfc |G7fc(i)/ifc(i)| 

Note that the constraint equation (Pf3l) contains higher derivative terms of Q like Q', 
which are not solved in our simulations. Hence it is impossible to evaluate the numerical 
accuracy of ( f43|) from our numerical simulation. In fact, if we rewrite the higher- 
derivative terms of Q by Q using Eq. (1581) . Eq. (|43|) becomes an identity. 

Fig. [5] shows that the relative error of Eq. ( 1421) is suppressed to less than 10~ 5 and 
Eqs. f )44|) and f|45|) are very precisely satisfied, indicating that our numerical scheme for 
solving the coupling system ( |25l) and ( |55l) is reliable. Note that the constraint equation 
( ]45i) has a different character from the others, which measures the numerical accuracy 
of the master equation (|25|) on the brane. 

5. Conclusion 



In this paper, we studied the evolution of the curvature perturbations in a brane-world 
inflation model where an inflaton is living on a single brane in a 5D AdS spacetime. 
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Figure 5. The numerical errors of constraint equations (|42|1 [upper-left] . 1|44|) [upper- 
right] and (|4"5)) [/ower] . These plots are results of calculations for k/a (ti)fi — 296. 

We used the Hawkins-Lidsey model of the inflaton potential which enables us to derive 
the background solutions analytically. We solved the full coupled evolution equations 
for the inflaton perturbations described by the Mukhanov- Sasaki variable and the bulk 
metric perturbations described by the master variable. This is the first numerical result 
for the evolution of scalar curvature perturbations in a brane-world that consistently 
takes into account the backreaction of the metric perturbations. 

We used the GN coordinate to describe the bulk spacetime that has a coordinate 
singularity. Then we are forced to introduce a regulator brane to cut the spacetime. 
We have checked that the evolution of curvature perturbations is insensitive to the 
location of the regulator. This is because the bulk metric perturbations decay fast away 
from the brane. We adopted the simplest possible initial condition for Q by taking 
Q = and Q = at an initial time. Then we fix Q so that the boundary condition is 
consistent. In order to check whether we can consistently neglect the coupling to bulk 
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metric perturbations described by J(O) in Mukhanov- Sasaki equation (1551) or not, we 
followed the subsequent evolution of curvature perturbations numerically. 

The evolution of curvature perturbations showed a suppression of the amplitude 
compared with the conventional 4D model even on sub-horizon scales. This means that 
it is impossible to neglect a coupling to gravity even on small scales due to the coupling 
to the higher- dimensional gravity through J(Q). The suppression of the amplitude may 
be understood as a loss of energy due to the excitation of the bulk metric perturbations 
as we took the initial condition that Q(y,ti) = and Cl(y,ti) = 0. On super-horizon 
scales, the curvature perturbations become constant, which confirms the fact that the 
constancy of the curvature perturbations is independent of gravitational theory for 
adiabatic perturbations [27] . 

Our result suggests that an usual assumption that the inflaton perturbations (the 
Mukhanov- Sasaki variable) approach to a free massless field on small scales cannot be 
applied in a brane world models on small scales at high energies. Then the spectrum 
becomes sensitive to initial conditions not only for the inflaton perturbations but also 
for bulk metric perturbations. In this sense, we should not take the result Fig. [3] at face 
value because this result is based on a special choice of initial conditions. In order to 
determine the initial conditions without ambiguity, it is required to quantise the coupled 
brane (inflaton)-bulk (metric perturbations) system in a consistent manner along the 
line of [28] . This is an outstanding open question to be addressed in order to derive the 
spectrum of the inflaton perturbations in a brane-world. 
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